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Abstract 

The physics of fc-core percolation pertains to those systems whose constituents require a minimum 
number of k connections to each other in order to participate in any clustering phenomenon. 
Examples of such a phenomenon range from orientational ordering in solid ortho-para H2 mixtures 
to the onset of rigidity in bar-joint networks to dynamical arrest in glass- forming liquids. Unlike 
ordinary {k = 1) and biconnected (k = 2) percolation, the mean field k > 3-core percolation 
transition is both continuous and discontinuous, i.e. there is a jump in the order parameter 
accompanied with a diverging length scale. To determine whether or not this hybrid transition 
survives in finite dimensions, we present al/d expansion for k-coie percolation on the d-dimensional 
hypercubic lattice. We show that to order the singularity in the order parameter and in the 
susceptibility occur at the same value of the occupation probability. This result suggests that the 
unusual hybrid nature of the mean field fc-core transition survives in high dimensions. 
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I. INTRODUCTION 



In a number of physical problems long-range order requires more than a single stranded 
path to propagate over long distances. One example of this is the propagation of quadrupo- 
lar order in solid (o — }l2)x{p — H2)i-a; (ortho-para H2) mixtures j^. This system can be 
reasonably modeled as a quenched, site-diluted lattice of electrostatic quadrupoles interact- 
ing via nearest-neighbor interactions on an fee lattice. Since the lowest state of two such 
quadrupoles at displacement r is one in which one molecule may assume any orientation in 
the plane perpendicular to r, it is clear that long-range order can not propagate down a long 
noncoUinear single stranded path. Therefore, to develop long-range order quadrupoles in an 
"infinite cluster" must have more than two nearest-neighboring quadrupoles. 
_Ariother such example of multi-path long-range ordering is rigidity percolation 

where each occupied site on a lattice has g degrees of freedom. The degrees of 
freedom of the site become fixed as more neighboring sites become occupied — one occupied 
neighbor constrains one degree of freedom. Therefore, in order to participate in the infinite 
rigid cluster, an occupied site must have at least g occupied neighbors. Here, again, is 
an example of a constraint on the minimum number of occupied neighbors giving rise to 
multiple long-range paths through the system. 

More recently, an analogy between multi-path percolation and the onset of elasticity 
in repulsive soft spheres as the packing fraction of the system is increased has been put 
forth by Schwarz et al. (SLC). In the analogy overlaps between particles correspond to 
occupied neighboring sites. To ensure local mechanical stability for each particle, d + 1 
occupied neighbors are required for each occupied site, otherwise the site is unstable and 
it is removed from the system (as opposed to running into other particles in the system). 
Here, the d + 1 constraint gives rise to multiply connected paths that eventually span the 
system. The onset of elasticity in the repulsive soft sphere system — a type of jamming 
transition called Point J — is thought to have implications for other phenomena such as the 
glass transition and the colloidal glass transition ^ 0] . 

The particular model of percolation called /c-core, or bootstrap, percolation 
turns out to be the relevant model of interest for such systems. For both the solid ortho-para 
H2 mixtures and the jamming system, fc-core percolation is an approximate description. 
However, for rigidity percolation, at least on the Bethe lattice, (/-rigidity percolation is 
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FIG. 1: Culling of a cluster for A; = 3 on a Bethe lattice with coordination number z = 4. The 
filled circles indicate branches which are k — 1 connected to infinity. Culling removes the bonds at 
a, b, and eventually at c, including the lower two sets of filled circles, but the rest of the cluster 
survives. 

equivalent to {g + l)-core percolation In fc-core percolation, each bond is independently 
occupied with probability p and vacant with probability 1—p. In addition, it is required that 
sites with less than k occupied neighbors should be made vacant. This "culling" operation 
proceeds recursively until all remaining unculled occupied sites have at least k neighbors, as 
illustrated in Fig. ^for a Bethe lattice where each site has z neighbors with z = 4. Such a 
model gives rise to "many" paths emanating to infinity from a single site for large enough 
p. 

Analysis of fc-core percolation on the Bethe lattice by SLC showed that the critical point 
has some unusual characteristics. As was known previously Q|, they found that the transi- 
tion at p = Pq — the p at which an infinite cluster appears — is a discontinuous one for k > 3, 
and the probability Pqo that a site be in the infinite cluster (after culling) exhibits the power 
law behavior: 

^oo = 0, p< po, 
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P + a{p-po) 



1/2 



P>PO- 



(1) 



More strikingly, SLC found a diverging correlation length and that an appropriate suscep- 
tibility has a power law divergence at the same threshold value of p at which the order 
parameter has a discontinuity. 

An immediate question then arises, is this unusual structure of the critical point charac- 
teristic of infinite d, or does it survive for finite rf? A clear answer to this question could be 
provided if a field theoretical formulation of the problem were available, but at present, no 
such field theory exists. (Although a field theory for the jamming transition based on force 
balance, as opposed to local mechanical stability, has been recently proposed [l^, it is not 
obvious that this model is in the same universality class as fc-core percolation in finite dimen- 
sions.) As for lattice models in low d spatial dimensions, until recently, numerical studies 
indicate that the fc-core transition is either continuous or pc = 1, i.e. nothing survives the 



culling process until the lattice is fully occupied [IJ, ll5[ . A case of the former is = 3 on 
the triangular lattice 1^, and for the latter, = 3 on the square lattice [IJ. In fact, for 
A; = 3 on the square lattice it has been proven by van Enter that Pc = 1 llSH. Moreover, 
it has recently been proven for hypercubic lattices that pc = 1 for k > d + 1 To date, 
the only numerical evidence for a k-core transition with a jump in the order parameter and 
yet a diverging correlation length is a directed k-core model with 10 nearest neighbors in two 
dimensions studied by SLC 20|. This result supports the notion that the unusual nature 
of the k-core transition found in the mean field survives in finite dimensions, though simple 
lattices, like the square lattice, are probably not the place to look for such a transition. 

Lacking a k-core field theory and the sparsity of finite d numerical results, we decided to 
implement an expansion in 1/d. As one sees for the Ising model, this expansion cannot be 
used to discuss critical exponents, because these exponents are independent of d for d > dc, 
where dc is the so-called upper critical dimension. For the Ising model, (ic = 4 j22, 23| and for 
ordinary percolation is = 6 • However, this expansion has been used to generate 

short series in l/d for the critical value of the coupling constant in problems such as the 
Ising model and self-avoiding walks, 3| for spin glasses, [^l for lattice animals, j3| and 
for ordinary percolation [2^. Thus, this expansion is ideally suited to answer the question 
or whether or not the unusual k-core mean field transition survives in finite dimensions. 
It should also be noted that the interpolation to continuous dimension implied by this 
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expansion is precisely the same j30| as that used by Wilson in his development of the 



renormalization group [2^. We will use this expansion to calculate corrections up to order 
for the critical coupling constant at which the discontinuity in the order parameter takes 
place and compare it to the critical value at which the susceptibility diverges. Since we find 
that these two threshold remain equal up to this order in 1/c?, we conclude that it is likely 
that this coincidence is robust and remains true at least for some range of high d. We 
note that these results are stronger than the results of Toninelli et al. [l9| showing that the 
unusual nature of the transition survives on Husimi trees (a Bethe lattice with finite loops), 
since, like the Bethe lattice, this structure has an infinite fractal dimension and therefore 
provides no information on the situation for finite d. 

Briefly this paper is organized as follows. In Section 11 we review known results for the 
Bethe lattice to fix our notation. Section III introduces the notion of perturbing the equation 
of state. Section IV introduces the concepts behind the 1/d expansion. Section V presents 
the perturbative 1/d corrections to the equation of state, while Section VI does the same for 
the susceptibility. In Sec. VII we summarize the conclusions which may be drawn from our 
result, that to order 1/d^ at large d, the k-core transition remains a hybrid transition and 
in Sec. VIII we discuss implications of our result for systems such as glass-forming liquids. 



II. BETHE LATTICE EQUATION OF STATE 

In this section we construct the self-consistent equation for the order parameter on the 
Bethe lattice for k-coie bond percolation. For this purpose we consider a rooted Bethe lattice, 
i.e. a lattice emanating from a seed site (as in Fig. 12)) in which the lattice is constructed by 
recursively adding a = z ~ 1 sites to each bond. To construct the self-consistent equation, 
the missing zth neighbor of the seed site is posited to survive culling. Therefore, the entire 
cluster will survive culling if each recursively added site has k — 1 outward bond connections 
to infinity. We therefore define the quantity P^o to be the probability that when we add a 
site to the cluster, that site is then has k — 1 outward bond connections to infinity which 
survive culling. The probability that some site on a Bethe lattice is in the /c-core can then 
be related to and has the same type of singular behavior as Poo for A; > 3 [3]. Accordingly, 
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we have the self-consistent equation of Poo as j?, 10 1 

m=k—2 



1- E 



which we write in terms of Q = pPoo as 

m=fe— 2 



[l-pPooT-'^ipPo, 



Q/p 



E 

m=0 



[a — m]\m\ 



-Q^{l-QY~'^ = ^{Q). 



Starting at p = 1, so that Q = 1, we consider the effect of reducing p and find that 



1 d$ 
p dQ 



As p is decreased, Q decreases until p reaches a critical value pc at which 
1 d^{Q) (a-fc + 1)! 



Q'-\l-Q) 



(2) 



(3) 



(4) 



(5) 



p dQ {a - k + l)\{k - ?>)\ 

Indeed, when the description of the transition is phrased in this way, it seems almost obvious 
that "turning on" finite spatial dimension will not invalidate Eq. (fT]). 



III. PERTURBATIVE CORRECTIONS TO ^{Q) 

Now we wish to perturb the equation of state (EOS). For instance, we may consider what 
terms appear in the expansion of the EOS which depend on the variables PaPbPcPdi where a, 
6, c, and d form a square. We will construct 

A<l> = $H - $5 , (6) 

where the subscript on $ indicates whether it is to be evaluated for the hypercubic (H) lattice 
or the Bethe (B) lattice. This expression includes the additional iterative term (which we 
only invoke once at order l/cr^ where a is 0{d)) for the H lattice and it also takes account 
of terms which appear on the B lattice but which do not have counterparts on the H lattice. 

To evaluate the derivative with respect to the p's in this expansion, we note that each 
occupied bond carries a factor of p and each unoccupied bond carries a factor of 1 — p. 
So for any bare diagram of h bonds, the derivative with respect to all its p's will involve 
a sum of the 2^ configurations of occupied and unoccupied bonds, in which an occupied 
bond carries a factor +1 and an unoccupied bond carries a factor —1. Thus, from a square 
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FIG. 2: Diagrammatic representation of the EOS. Left: the unperturbed EOS based on a seed site, 
S, with a emerging bonds. Right: the perturbation to the EOS from the insertion of a square, so 
that now the seed site has a — 2 single bonds and an attached square. 




FIG. 3: Contributions from the 16 subdiagrams (of diagram a) which result from either occupying 
(full line) or not occupying (dashed line) each bond of the square on the H lattice. Greek letters 
label sets of subdiagrams which cancel one another. 

(which is the configuration giving the leading correction at relative order 1/cr^) we generate 
16 subdiagrams (the first of which is just the one for which all bonds are occupied and the 
last of which is the one for which all bonds are unoccupied). Accordingly, the derivative 
with respect to the four p's which form a square (the leading correction to the B lattice) 
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FIG. 4: As above, but for the B lattice. Vertices which are cfose to one another would coincide on 
the H lattice. 

gives rise to Figs. 1 for the H lattice and to Figs. 2-5 for the B lattice. 

If f{x) denotes the contribution to A$ from square insertions, which we call $sq) then 

$sq = f{a) - 2f{b) - 2f{d) - 2/(0 + 2/(j) + 2f{k) - f{m) + 2f{n) - 2f{q) + 2/(r) (7) 



where we noted that /(/) = f{h) = /(/) and f{g) = f{p). This equation is represented in 
Fig. [7| The factors of 2 take account of the fact that some topologies occur in two equivalent 
realizations. When we calculate the contributions / we must include not only the factor 
PaPbPcPd = but also the number of ways (which depends on cr) the bonds making up the 
square can be selected out of the a available bonds. Since this factor is of order o"^, we see 
that $sq is of order p^a"^ which we will see is of order a^^. As will also be seen later, if we 
were only evaluating these corrections to leading order in 1/a, then we would not need to 
differentiate between diagrams (b), (j), and (r), or between (d), (k), and (n). Accordingly, 
at order 1/cr^, when we consider the analogous contribution from hexagons, ^Hex, we do not 
have to be explicit about the configuration of the vacant bond(s) and we therefore have the 
hexagon insertions shown in Fig. |H| 
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FIG. 5: As above, but for the B lattice. 




FIG. 6: As above, but for the B lattice. 
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FIG. 7: Diagrammatic representation of Eq. [7J Note that there are no diagrams with two unoccu- 
pied bonds. This is because for every H lattice diagram there is a B lattice equivalent with a "cut" 
between the two unoccupied bonds resulting in a cancellation. 



A. Percolation 

We can apply the above formulation to calculate the corrections to the critical concen- 
tration pc in the ordinary bond percolation problem, and check the results by comparing 



to those obtained by Gaunt and Ruskin [29|. For percolation, a diagram contributes to the 
probability of the seed being connected to infinity if one or more of its vertices is connected 
to infinity. So, for the rooted Bethe lattice the EOS at the transition is 

P = 1 - {1-pPy . (8) 

(The only nonpercolating case is if all vertices are not connected to infinity.) Since the 
percolation transition is continuous, the critical concentration is found by expanding the 
EOS in powers of P. At linear order the solution requires that 

1 = ap,. (9) 

For hypercubic lattices the contribution to •I'sq from Eq. ((Zj) is written as 



^^^ = ^{cy-\Yp 



[1 - (1 - ppy-'^ii - p"-')^] - 2[i - (1 - pp)""'(i - p'^-^f] 

-2[1 - (1 - pP)"-2(l _ p-l)3] _ 2[1 - (1 - P)2(l _ pPy-'^il - P--1)2] 
+2[1 - (1 - P)(l - pP)"-2(l _ p-l)2] + 2[1 - (1 - P)(l - pP)"-2(l _ p--l)2] 
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(a) 



-2 



(b) 



-2( 



(c) 



-2 



(d) 



-1 




(e) 



+2 



(f) 



+2 




(g) 



+2 




(h) 



FIG. 8: As in Fig. [3 contributions to ^{Q) from hexagons on either a H or a B lattice. Diagrams 
(f), (g) and (h) correspond to sums of graphs on both the H and B lattice. For example, the square 
equivalent to hexagon diagram (f) is the sum of diagrams (b), (r), and (j) in Fig. [7| Diagrams 
similar to these give the contribution to the susceptibility, Th,a, due to hexagon insertions as we 
will eventually show. 



-[1 - (1 - Py{l - pPY~\l - P' 



2[1 - (1 - P)(l - pPY~\l - P 



1(T-l\2l 



-2[l-(l-P)(l-pP) 



cr-l. 



l-P 



cr-l\3l 



2[l-(l-pP)'""^(l-P'' 



(10) 



where the terms between each small square bracket denote diagrams (a), (b), (d), etc. from 
Fig. 13 and 



P 



cr-l 



l-{l-pPY~^ + 0{a~'). 



(11) 



P"^"^ is the probability that a site has at least k — 1 connections to infinity out of a — 1 
bonds attached to the site. A generalization of this quantity will be introduced later on. 
Then $sq can be simplified to read 



1 



<^'Sq=2(^-l)V 



2P(1 - pP)'"-^(l - P^-^f + 3(1 - pP)'""^(l - P^-i)3 



+3(1 -pP) 



a-2i 



P 



(7-l\2, 



l-Pf - 6(1 -pP) 



<T-2l 



:><T-1\2, 



l-P) 



:i2) 



The factors of (1 — pPy and (1 — pPy ensure that the bonds attached to the seed 
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FIG. 9: Attaching a square to a site. 

site but not in the square are isolated from infinity, and the factor A'^sq = (o" — 1)^/2 counts 
the number of possible squares. We construct N^q using Fig. IHl as follows. Say the bond to 
which we wish to attach the square is in the xi direction. Then there are two cases: either 
the square involves a bond in xi direction or it doesn't. In the first case the square involves 
the Xi bond and one of 2c/ — 2 = a — 1 choices for the other axis of the square. The xi bond 
can either be the first or the second leg of the square, but we ignore this degeneracy because 
it simply corresponds to traversing the square in different senses. So for this case there are 
(cr — 1) configurations. In the second case the square involves two out of — 1 dimensions, 
so there are 4((i — l){d — 2)/2 = (a — 1)((t — 3)/2 such configurations. In all 

iVsq = (a - 1)2/2 . (13) 

This number would be twice as large if the sense in which the square was traversed mattered. 
For hexagons the calculation is much simpler because we do not need to keep track of the 
configuration of vacant bonds, which forces us to distinguish between P and P'^~^ and which 
leads to corrections of relative order 1/a. Using Fig. |H1 we write the contribution to A$ 
from hexagons in leading order of 1/cr as 



5(1 - P)^ (1 - [1 - P]) + 2(1 - pP){l - Pf (1 - [1 - P]) 



= -2ay (1 - pPy-'P{l - P)^[5 + 2(1 - pP)] . (14) 

Here the factor 2a^ is (to leading order in l/o") the number of ways of attaching a hexagon 
to a site. Note that p ^ 1/cr + 0{a~^) so that to the order we need, P'^^^ ~ P[l — (l/a)]. 
Also note that for this continuous transition we only need the contribution to A$ which are 
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linear in P. Thus we set A$Hex = — 14P(T ^ and then 

A<l> = 



p r. „ P 



1 - 



[2 + 3[l-(l/a)]-14-. (15) 



2^2 

Finally, using 1 = apc + dA^/dP, we get 

ap, = 1 + (5/2)(T-2 + (15/2)(T-3 + Oia-^) , (16) 
in agreement with Ref. 0. 

IV. EXPANSION IN POWERS OF 1/a 

Before getting into the calculation we should indicate how the various variables are to be 
expanded in powers of l/a. Consider Eqs. ^ and (0). Set 

Po = E /^"^"" ' = E ' (17) 

n n 

where the subscripts "0" indicate that the quantities are those of the Bethe lattice solution. 
Although we will not invoke the actual values of the coefficients in these expansions, we will, 
for illustrative purposes, determine the leading coefficients, ai = a and /?i = P, where 

m=k-2 m 

and 

1 = /3a'=-2e-" . (19) 
By eliminating /3, we obtain an equation which determines a: 



m=k—2 



+ E ^- (20) 



m=0 ^ 



For k = 3, numerical evaluation yields the approximate values a = 1.793282133 and (3 = 
3.350918872. As we shall see, it is more convenient to express results in terms of po and Qo 
because our main aim is not to explicitly obtain an expansion for these quantities in powers 
of 1/cr, but rather to determine whether or not the singularity in the EOS state coincides 
with the singularity in the susceptibility. 

Now we see how the critical coupling constants Pc{d),Qc{d) for hypercubic lattices are 
obtained after A$ has been evaluated. We write 
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R 



a- 1 



k-2 



R 



t 

a-2 
k-2 



SEED 




FIG. 10: Definition of i?^ (in panel a) and P"" ^ (in panel b). Each site has cr + 1 bonds, 
and set 

p = Po + Ap, Q = Qo + AQ . (22) 

In this analysis we note that A$ ~ a^^, so that to obtain results to order we need only 
consider terms linear in A$, Ap, or AQ. Then we have 



Qo ^ AQ QpAp 
Po Po Po 



which gives 



Po 



Also 



Ap=-^A^Qo,po) . 
vo 



]__Ap ^ rf$o(Qo)_^^g^ , rfA<l>(go,Po) 
Po Po dQ dQ'^ 



which gives 



AQ 



dQ 






-1 - 


_dQ\ 




'd^% 


-1 r 


_dQ\ 





Po dQ 



dQ 



A$(go,Po) dA^ 



Qo 



dQ 



(23) 



(24) 



(25) 



(26) 



We introduce the following Bethe lattice quantities. First, i?^"*" is defined to be the 
probability that out of o" — m available bonds, exactly r are — 1 connected outward (i. e. 
away from the origin or seed of the cluster) to infinity and we set R^^^ = Rr- (It is also 
convenient to set Rk-2 = R)- This definition of R is illustrated in Fig. ^| Similarly we 
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define P'^-"^ to be the probabifity that a site is at least k — 1 connected to infinity through 
the set of (T — m of bonds and P"^ = P. Since Q is of order l/cr, we see that P'^~'^ is of 
order (assuming m is of order unity). One can hkewise show that under this same 

assumption R^^"^ is also of order unity. Note when the superscript on P is of the form a — n 
then the superscript indicates the number of bonds, as defined above. When the superscript 
is is a purely numerical value like "2", then this indicates an exponent: P^ — P x P. 

We now express all the above Bethe lattice quantities in terms of the canonical variables, 
R, P, and Qq = p^P. We have 



m=k—2 



pa- 



1- E 



(a-r)! 



™=o (o--r-m)!m! 
from which we obtain, correct to first order in l/cr that 



(27) 



(28) 



Also, 



R. 



Qo 



Rr, 



{a-l-r)\ 



m\{a — 1 — r — m)! 



:i - Qo 



\a—l—r—m 



l+rQo- 



rm 



(29) 
(30) 



and 



dP 



(a-l)! 



(m- l)!((7-m- 1)! 



Qr'(l-Qo. 



\u—m—l 



(a-l)! 



(31) 



m!((7 — m — 2)! 



Q^(i - Qo) 



u—m—2 



1 



(a-l)! 



(m-l)!((7-m-l)! 

ia — m)Rm-i — (a — m — 1)R 



{a-iy. 



m\(a — m — 2)! 



Wi-Qo) 



a—m—l 



i-Qo 

(a-m + Q) {Rm-i - Rm) + Rm + 0{l/a) . 



(32) 



V. 1/a EXPANSION FOR THE EQUATION OF STATE 



In this section we will implement the 1/a expansion for the EOS. Note that we are 
considering the effect of having a square (or hexagon) of bonds with one vertex at the seed 
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FIG. 11: Configurations of diagram (a) of Fig. IS] which have two paths to infinity. The heavy line 
is a bond which is A; — 1-connected to infinity and carries a factor P^^^. The dashed hne carries a 
factor Rk-2 which may survive culhng if it is connected to two hve bonds. 

site. (Fig. El ) All structures emanating from this square (or hexagon) or from the other a— 2 
bonds which intersect the seed site, may be assumed to be tree-like because the occurrence of 
more than one loop only influences terms of relative order and we do not consider this 
order. Accordingly we implement Eq. ((Tj). For squares we have f{n) = [(cr — l)^p^/2]5/(n), 
where n indicates the diagram as labeled in Eq. ((7j). Since the insertion of the square (or 
hexagon) can contribute a maximum of two more paths to the boundary (or infinity) at the 
seed site, we will break up diagram (a) of Fig. 3 into factors associated with having two, 
one, and zero paths to infinity as illustrated in Figs. [TTJ [T21 andHSl respectively. Thereby 
we find 



Sf{a) = (rIzI + RUl + P"'^) [{P^'^f + 2i?(P"~i)2 + m^P''-^ + E? 



+ [Rk-l + P""^][2 + 4i?]P"-i[l - P"-i - R] 



+ P 



cr-2 



1 - P''-^ - Ry + 2P(1 - P''-^ - Ry + 3P^(l - P"-^ - R) 



(^pa-iy ^ 2R{P''-y + SR^P"-^ + R 
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FIG. 12: As in Fig. 111! configurations of diagram (a) of Fig. |31which have one path to infinity. The 
symbol K denotes a bond which is definitely culled and therefore carries a factor (1 — P°"^^ — Rk-2)- 

.K /v . K 









FIG. 13: As in Figs. 1111 and fHH configurations of diagram (a) of Fig. |31 which have no paths to 
infinity. 



P0--2 



1 _ _ jiy + 2R{1 - P 



0--1 



Rf 



+3R\l - P"-^ -R) 



(33) 



In constructing this expression we noted that if the square had two paths to infinity, then 
the remaining cr — 2 bonds emanating from the seed site had to have at least k — 3 paths 
to infinity which we take into account with the factor of the form Rks + Rk-2 + P- {P is 
the probability of having at least k — 1 paths to infinity.) The factor Rk-2 + P and P take 
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account of having at least A; — 2 or A; — 1 paths to infinity, respectively. Similarly, for the 
other diagrams we have 

Sf{b) = [P"-'' + irkZi]P''~'{l + R + R'') + P''~''[l- p^-'ii + R + R")] 

= P"-^ + RlZlP''~\l + R + R^) (34) 



6f{d) = [P^-' + RlZl + RlZl]{P''-'n^ + R) 

+P'^-2[1 - P'^-l][l - P'^-l(l + P)] 

= P'^-^ + P^^' + RlzliP'-'^il + R) 
-RlZl[l-P''-'][l-P'"\l + R)] 



(35) 



Sf{i) = P'^-^ll - P][l - pr'(i + R)- PR 



[1-P][P'^-^(1 + P) + PP2] 



+P[l-Pr'(l + i?)-^i?')] 
+(P^I^ + P^Z^ + P'^-2)P[P^-^(1 + P) + pp2] 

+ P"-^ + P^^p^pa-l^^ + P) + PP2] 



po--2 



RU{i-P) 



l-P^-\l + R)-PR' 



(36) 



5/(j) = [P— 2 + P-I2][p— + P) + PP2] + P''-2[l - P''-^(l + P) - PP2] 

= P'^-^ + P^i^jP'^-i^^ + P) + Pp2] (37) 



5f(k) = [P'^-2 + Pri + P^Z3^]PP'^-^(l + P) 



+[P''-2 + p; 



.(T-2] 



p— 1[1 _ p](l + p) + p[l _ P— ^(1 + P) 



+P''-'[1 - P][l - P''-'(l + P)] 



= P"-' + RlZi + RlZiPP"' {I + i?) 
-P^Z2[1-P][1-P'^-^(1 + P)] . 



(38) 
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Sf{m) 



a-1 



+2[P 

pa-2 _^ 



R 

"-"^ + R 



(T-2 

k-2 



RPf 

][1 - P"-^ 



RP][P^-' + RP] 



a-2 ^p<T-2j^pa-l _^^p|2 



fe-2 



RPf 



(39) 



6 fin) = P^-^ + Rlzl + RlzlP''-'[P''-' + RP] 
-Rlzlll - P""^][l - P""^ - PR] 



(40) 



5/(g) = P 



cr-l 



I- P''-\l + R + R') - R'P 



+{R + P 



P''-\l + R) + R') + R-'P 



P"-^ + R 



P'^^Vl + P + P^) + P^P 



5 fir) = [P"-i + P]P"-i[l + P + p2]+P'^-i[l-P'^-i(l + P + p2)] 
= P"-^ + PP'""^ [1 + P + p2] . 

With these results we are now ready to evaluate Eq. ((7j) and it reads 

ia - 1) V '' 



Sq 



RlZi[R + 3R^P''-^ - 3P^P^] - 2P^P 



+RlZl[R^ - SR'P^-' + 3P"P 



>2 d2i 



(41) 



(42) 



(43) 



Since the contributions from hexagons only need to be evaluated to leading order in l/cr, 
the configuration of vacant bonds is irrelevant and we may omit the superscripts on P and 
P. Accordingly, for hexagons it is convenient to classify (as indicated by the subscript) the 
value of m, the number of paths to infinity. For diagram x we write 



fmix) = 2a'' 5 f mix 



(44) 



where 2a^ is the number of hexagons (to lowest order in l/o") that can be attached to a 
bond and 



5/2(0) = p2[l + 2P + 3P2 + 4P^] + 5P^P + P^ 



(45) 
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Sfi (a) = 2P[1- P - R][1 + 2R + 3R^ + 47?^] (46) 

Sfo{a) = [1-P-Rf[l + 2R + 3R^ + 4R^]+5R'^[1-P-R] (47) 

Sf^{b) = P[l + R + R^ + R^ + R"^ + R^] (48) 

5fo(b) = [1-P-R][1 + R + R^ + R^ + R'^]+R^[1-P] (49) 

(5/2 (c) = P2[l + + + + i?^] (50) 

5/i(c) = [1 - P]P[1 + + i?' + i?' + i?^] 

+P[1 - P - P] [1 + P + p2 + P^] + P[l - P]P^ (51) 

5/o(c) = [1-P][l-P-P][l + P + P' + P3] + [1-P]2p^ (52) 

5f2 (d) = P' [1 + P + P' + P^'] [1 + P] (53) 

Sf^{d) = P[l- P - RP][1 + R + R^ + R^] 

+P[1 + P] [(1 - P - P) (1 + P + p2) + (1 - P)p3] (54) 

5fo{d) = [1-P-PP][(1-P-P)(1 + P + P') + (1-P)P'] (55) 

Sf2{e) = [P(l + P + P2)]2 (56) 

5/i(e) = 2P[l + P + p2][(l-P-P)(l + P) + p2(l-P)] (57) 

6fo{e) = [(1 - P - P)(l + P) + p2(i - P)]2 (58) 

(/) = P[l + + i?2 ^ i?3 ^ (5g^ 

5fo{f) = [(1-P-P)(1 + P + P' + P3) + (1-P)P^] (60) 
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Sf2ig) = P^[l + R + R' + R^] 



(61) 



6h{g) = P{l-P)[l + R + R'' + R'] 

+P[{1~P-R){1 + R + R^) + {1-P)R^] (62) 

Sfoig) = {1 - Pm + R + R'){1 - P - R) + {1 - P)R'] (63) 

Sf2{h) = P\l + R + R^){l + R) (64) 

Sf^{h) = {1-P-RP)P{1 + R + R^) 

+P{1 + R){1- P - RP - R^P) (65) 

5/o(/i) = (1 - P - RP){1 -P-RP- R^p) , (66) 

where the numbering of contributions is as in Fig. |H1 Thus for the contribution of hexagons 
to the EOS (indicated by the superscript "H"), the sum over diagrams gives 

$f = Ua) - 2/2(c) - 2h{d) - h{e) + 2h{g) + 2h{h) 

= 2aY (^-5R^Pl, + 5R*Poo + R'^ • (67) 

$f = A (a) - 2/i(6) - 2/i(c) - 2A(d) - A(e) + 2A(/) + 2/1(^7) + 2/i(/.) 

= 2aV (^WR^P^ - (lOi?^ + 2i?^)Poo j , (68) 

and 

< = Ma) - 2/0(6) - 2/o(c) - 2/o(d) - /o(e) + 2/o(/) + 2/0(^7) + 2/o(/i) 

= 2aY (^-5R^Pl + (5P^ + 2p5)Poo - i?'^ • (69) 

It is a check on our resuhs that J2n — 0- Now we match each of these contributions to 
possible configurations of the other a — 2 bonds from the seed site. Thus, from hexagons we 
get 

$Hex = P<^0 +iP + R)^? + {P + R + i?fc-3)$f 

= Pfc_3$f - R<^^ . (70) 
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VI. 1/ct expansion for the susceptibility 



A. Formulation of the Susceptibility 

In this section we consider the 1/a expansion for the two-point susceptibihty, Xij- 
For fc-core percolation we calculate this quantity in the ordered phase, because for the Bethe 



where Ui is an indicator variable which is unity if the site i is in the fc-core and is zero 
otherwise. Also the angle brackets indicate an average of configurations of bonds in which 
each bond is independently present with probability p and absent with probability 1 — p. 

Now we formulate the 1/a expansion for the ordered phase of k-core percolation. On the 
Bethe lattice the long-range part of this correlation function comes from the probability of 
configurations in which two sites are only in the k-core by virtue of the presence of a path of 
occupied bonds connecting sites i and j. (Of course each site on this path must belong to the 
fc-core.) In this connection, it is important that the sites not be in the fc-core if any bond is 
removed. Such contributions are either already counted in lower order or are canceled when 
the second term in Eq. ()71|) is subtracted off. Thus the configuration in panel (a) of Fig. 
is an allowed configuration contributing to Xijj but that in panel (b) can not be extended. 
In Ref. 3) this consideration led to equivalently restricting the sum to the "corona." Thus, 
as the path is progressively lengthened, if we reach a point where the origin of the path 
is certainly in the k-core no matter how the path is extended, then this path is said to be 
"truncated," and is discarded as not contributing to the singularity in x- 

For the present work we obtain the results of Ref. Q as follows. For a path to satisfy the 
"corona" constraint, it must consist of a path for which each vertex is — 2 connected to 
infinity if the bonds along the path are not considered. With this construction one sees that 
at each vertex, whether the whole structure is or is not ^-connected to infinity depends on 
what happens further down the path. Thus, for the Bethe lattice each vertex in the path 
leading from site i to site j carries a factor pR'^Zl = pRk-2 = pR- When Xij is summed over 
j, one obtains a geometric series in the ratio, r, given by 



lattice there are no finite A;-core clusters. Here we define x = Xij^ where 




(71) 



r = 



(^pRk-2 , 



(72) 
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FIG. 14: Configurations for k = 3 illustrating "truncation" of diagrams. Here the side branches 
are assumed to be /c — 1-connected to infinity. In the top panel a configuration is shown where the 
k-cove depends on the state of the further bonds indicated by filled circles. In the lower panel, the 
configuration is /c = 3 connected no matter what happens further down the chain. Since the fc-core 
does not depend on further bond occupation probabilities this diagram can not be extended. 

which coincides with the singularity in the EOS 

A convenient way to describe the 1/a expansion as applied to the susceptibility, is to 
regard the 1/a corrections as a vertex correction, as illustrated in Fig. El for insertions of 
a square. This vertex correction amounts to replacing Rk-2 by Rk-2 + Ai?, where AR is 
the contribution from the dashed box (which replaces each vertex). Analogous corrections 
arise at relative order from insertions of hexagons. However, insertions involving two 
or more loops lead to corrections of order l/a^ and higher and are beyond the scope of 
the present calculation. (Note that for the last diagram, SQ3, we need an extra factor of 2 
because the two sides of the square are inequivalent, so that it occurs (cr — 1)^ ways.) The 
diagrams for a square insertion into a chain diagram for the susceptibility are shown more 
explicitly below. 
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FIG. 15: The factors associated with the dashed box in the first three diagrams renormahze the 
vertex (the filled circle) in the last diagram. We label these contributions SQl, SQ2, and SQ3, 
respectively. 



The condition for a divergent susceptibility is that the perturbed ratio r be unity, or 



r = apR + apAR 



(73) 



We now evaluate the ratio r at the singularity in the EOS, so that p and Q are given by 



Eqs. 



and 



Working to leading order we write 



r = cr(po + Ap) 



R{Qo) + ^AQ 



+ apAR 



^ Ap 1 dR 
Po R dQo 



Qo R dQo 



a 



dR 
dQo 



A<^ dA<^ 



Qo dQ 



apAR 



1 + 



AR 



adQ 



R-' 



(74) 



Here we noted that for the Bethe lattice apoR{Qo) = 1 and d^/dQ = aR. There are 
now three scenarios, depending on whether the ratio r is a) less than, b) equal to, or 
c) greater than unity. In case a) the discontinuity in the EOS preempts the divergence 
of the susceptibility and the transition is a conventional first order transition, except for 
a fractional power law in the EOS for p above threshold. Case b) (consistent with the 
expansion up to order 1/cr^) indicates that the coincidence of the singularities in the EOS 
and the susceptibility is robust and survives for large but finite spatial dimension. 



24 




FIG. 16: Diagrams for Square #1 on cither H or B lattices. Diagram (k) occurs on the Bcthe 
lattice but must be subtracted off because its counterpart on the hypercubic lattice has a different 
topology. 

B. Diagrams for the Susceptibility 

Now we consider the effect of inserting squares or hexagons into the Bethe lattice diagrams 
for the susceptibility. For squares we write 



Arsq = (a - l)2^Ai?/i?fe_2 , (75) 



where the square insertions are discussed below. Similarly, for hexagons we write 

ArHex = 2pV3Ai?/i?fe_2 , (76) 
where the hexagon insertions are also discussed below. 
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_ 4 f (b) > _ 2 i( (c) 



-4 




+ 4 < (e) / +4 




FIG. 17: Diagrams for Square #2 on either H or B lattices. 



(a) 



(b) 



+ 2 




-2 /(d) 



- 2 



(e) 



X- V 



-2 



(f) 



-2 



(g) 



+ 2 



(h) 



+ 2 



(i) 



G) 



+ 2 



(k) 



FIG. 18: Diagrams for Square #3 on either H or B lattices. 
C. Percolation 

As an example we first carry out the calculation for the two-point susceptibility of perco- 
lation and we confine our attention to the disordered phase, i.e. Poo — 0. (The calculation 
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for the ordered phase is quite similar.) First, consider the squares. The sum of the con- 
tributions from diagrams (a) through (j) of Fig. [TBI vanishes. For diagram (k) there are 
^Sq = ~ 1)^/2 ways to attach the pseudosquare (which can be traversed in two senses). 
So the dashed box of diagram (k) carries the factor 

Ar#i = -pV^^^l- (2/a)] . (77) 

The contribution from Fig. El is Ar*^ = —{a — l)^/cx^ and that from Fig. ^lis Ar*^ = 
— (1/2) (a — 1)'^/(T^. So the total contribution from squares is 

Arsq = -(5/2)a-2 + (l3/2)a-^ (78) 

The contribution from hexagons is 

ArHex = -14a-3 , (79) 

When these contributions are summed, one finds that the threshold for the divergence of 
the susceptibility agrees to order with that in Eq. (fT^ found from the EOS and 

with the previous work of Gaunt and Ruskin 2^. Of course, this was hardly surprising, 
but the present calculation illustrates the more complicated calculation needed for k-core 
percolation. 



D. fc-Core Percolation, Squares 



We now consider the contributions to the susceptibility from insertions of squares as in 
Fig. El We calculate the factor AR associated with the dashed box. In the results given 
below we include the prefactors (e. g. ±2 or ±1) written in Fig. El We find 



AR{a) 



^(^ - 2) [ 



2a 



R 



fc-4 



2R{P- 



(T-l\2 



+ 3R'P 



2 ncT-l 



+R 



cr — 3 

fc-3 



2P'^-i(l _ p'^-i -R)+ ARP^-^il - P 



R' 



R) 



1 _ p'^-i _ p)^(i + 2R) + m\i - P"-^ - R) 



(80) 



where x = (cr — l)^p^ and we used A^^sq = (a — 1)^/2. Also, the factor {a — 2) /a takes account 
of the fact that the line leaving the square has only a — 2 choices because it has to avoid the 
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square. This factor appears for the other diagrams, except that for (b) and (f ) the counting 
is less obvious and is discussed in Appendix A. We have 

^RlZliP^-^l + R + R^) + R'^P] \ , (81) 



AR{c) = -x^^^l Rlzlil - P)[l- P''-\l + R)- R^P] 



+RlZlP[l - + R)- R^P] 

+Rlzl{l - P)[P"-'(1 + R)+ R^P] 



+RlZiP[P''-\l + R) + R'P] \ , (82) 



AR{d) = -x^^^-^IrIzU^ - P"-' - RP^ 



2a 

+2Rlzl[l - P"~^ - RP] [P""^ + RP] + PfeZ4 [P""^ + PP]^ 1 , (83) 



AP(e) = -x ^"" '^H rIzI[1-P''-\1 + R + R')]+RIzIP''-^[1 + R + R']] (84) 



AP(/) = x^^-^IrIzI[1 - P^-^l + R + P')] + Pr3^""'[l + R + R'']] (85) 



. .^{naii - P-d + R) - R^P) . ^aiP-d . R) . R^R]} (86) 

AR{h) = -x^^^|pr|(P'^-^)'[l + P]+Pr3^[l-P^-^]P^-^[l + P] 

+P^Z^P'"-^[1 - P"-^ - RP"~^] 

+RlZl[l - P"-^] [1 - P^~' - RP"-^] \ , (87) 



+P^Z^P''-^ [1 - P'^-l - RP] + P^I^ _ pa-l| _ pa-1 _ i^p] I , (88) 
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AR{j) = x^^-^i^RlzlPP''-^[l + R] + Rrl[^-P]P''-^[l + R] 

+RlZlP[l - P"'^ - RP"'^] + Rlzl [1 -P][l- P""^ - RP"'^] I . (89) 

Also 

Ai?(fc) = -i?y (a - 1)2 . (90) 

So, in all from Fig. El we get 

Sj=,Ai?(y) = ^^!^^PR\RIZI - Rlzl) + ^(^- 2)^r4 ^_3^2p2 ^ 3^2p<x-l ^ ^3] _ ^^5 



+ ^^^^i?r|[6i?'i^' - 6i?2p^-i] + '^^^^^—^Rlzl[m\P''-^ - P^) - i?tpl) 

2(7 2(T 

We next consider the contributions from the diagrams of Fig. El We get 



Ai?(a) = ^i^_il|2i?(i?^Z22)2[l-P'^-i-i?] + 2i?P'^-i(i?r3)' 

+p2(pr3)' + R'RT-lRrS^ . (92) 

In Fig. El we show the evaluation of the last term in this result. Note that for a diagram 
with a loop, it can matter which way the diagram is entered. 



AP(6) = -2^— ^xp2 (rIzI[P''-^ + RP] + Rlzl[l - P"'^ - RP] \ , (93) 



AP(c) = -x^-^^RlRlzlP + RlZl{l - P)f , (94) 



AR{d) = -2^^^^^RRlzl[RlZlP-' + RlZlil - P-')] . (95) 



AR{e) = 2x^R'[RlzlP^-' + ^^'(1 - P"'')] ■ (96) 



ARif) = 2^^^^^RRlzl[Rr-lP + Rr-l{^ - P)] ■ (97) 
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FIG. 19: Diagrams for Eq. ()92() . In the left-hand diagram the bottom site will be culled if the 
bond "x" is culled. In that case the entire cluster will be culled. So the stability of this cluster 
depends on the state of bonds below "x," and therefore this diagram can be extended. In contrast, 
the right-hand cluster survives culling no matter what may be the state of bond "x." Hence this 
diagram truncates and we discard it. 

So, in all we get from the diagrams of Fig. [T7| 

a a 

+^-^^:^R'[Rr-l - Rr-imzl + ^rt-I] . (98) 

For the diagrams of Fig. we have 

+2RlzlRlZlP'''^{l - P"'^ - R) 

HRl-lfi^ - P"'' - - P""' + R)Y (99) 

as is illustrated in Fig. 1201 



For the other diagrams in Fig. we find 

Ai?(6) = -x ^''~^\ RlzlfR^ , (100) 



a 



AR{c) = +2xR-'Rlzi . (101) 



AR{d) = -2xR'[Rlzl{l - P) + RlZlP] , 
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(102) 




FIG. 20: Evaluation of Ai?(a) of Eq. 



(103) 



Ai?(/) 



-2x- 



a-1 



a 



Rr_l{l - P) + RlZlP 



j^CT— 2ri DO"— 1 



[1 - P'^"^ - RP] 



+Rl-4[P"-' + RP] 



(104) 



Ai?(<?) = -2x^^^Rlzl 



a 



(105) 



AP(/i) = 2x 



a- 1 



cr 



.0--2 
'A:-2 



P^Zall - - RP] + ^rs l^"' + RP] 



(106) 
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AR{i) = 2xR 



RlZiil - P"'' - RP''-^] + RlZiil + R]P 



-X 



''^^[Rt-l{l - P-') + RlzlP-'r 



(107) 



(108) 



AR{k) = 2x^[Rlzl{l-P 



o■-l^ 



RlZiP^-'][Rlz',{l -P)+ RlZiP] ■ (109) 



So, in all from the diagrams of Fig. ^Jwe get 

a - 1 



S^^,Ai?(2/) = AxR'P[Rr4-Rr4]+x 



a 

cr-l r>2 



{RlZiY[-2R' + 2R{P''-^ - P')] 



+ {RlZlf[R^ + 2i?(P'^-l - P2)] + ^^-2^^-2[^2 _ 4^(p<x-l _ p2)| k 



110) 



E. A;-Core Percolation, Hexagons 

Now consider the diagrams of Fig. |H1 but for the susceptibility instead. Here and below 
the vertices (or vertex) which are part of the chain are indicated by filled circles. With 
y = 2a^p^ we have (with R = Rk-2) 

AR{a) = y\ Rk-i[P\l + 2R + 3R^ + AR^) + 5R^P + R^] 



+2Rk-3P 



1 - R^ - P{1 + R + R^ + R^) + R{1 - R^) - RP{1 + R + R^) 



+R\1 - R^) - R^P{1 + R) + R%1 -R)- i?^P 



+R 



(1 - P - P)2(l + 2P + + AR^) + 5P^(1 -P-R) 



(111) 



AR{b) = -2y{Rk-3P{l + R + R^ + R^ + R'^ + R^) 



+R[1 - P{1 + R + R^ + R^ + R^ + R^)] 



(112) 



AP(c) = -2y|Pfe_4P^(l + P + P^ + P^ + P^) + Pfc-3^(1 - P)(l + P + P^ + P^ + P^) 

+Rk-3P[l - P(l + P + P' + P' + P')] + 

P(l-P)[l-P(l + P + p2 + P^ + P^)]l (113) 
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AR{d) = -2y|i?fc_4P'(l + R){1 + R + R'' + R') + Rk-3P{1 + - P(l + + i?' + R^)] 
+Rk-3P[l - P{1 + R)]{l + R + R^ + R^) 

+R[1 - P(l + R)] [1-P{1 + R + R^ + P=^)]| (114) 
AR{e) = -y\Rk-4P'^il + R + RY + '^Rk-3Pil + R + R^)[l- Pil + R + R^)] 



2m2 



+R[1-P{1 + R + R')] 



(115) 



AR{f) = 2y{Rk-3P{l + R + R^ + R^ + R^)] 



+R[1 - P{1 + R + R^ + R^ + R^)] 



(116) 



AR{g) = 2y^Rk-4P'il + R + R' + R-') + Rk^^Pil - P){1 + R + R' + R') 
+Rk~3P[l -Pil + R + R^ + R:')] 
+R{1 - P)[l - P(l + i? + + R^)] I 



(117) 



AR{h) = 22/|Pfc_4P'(l + P)(l + R + R^) + Pfc-3P[1 - P(l + P)](l + P + P' 
+Pfc_3P(l + P)[l - P(l + P + P')] 
+P[1-P(1 + P)][1~P(1 + P + P')]| . 

In addition we have the analog of diagram (k) of Fig. ^1 

AR{i) = -2yR^ , 

We call the sum of these contributions Tha- We have 



(118) 



(119) 



Pfc-4[-5P'P^ + 5PP^ + P""] + 2Pfc„3P[5P^P - 5P^ - P^] - 5P^p2 



+ (5P^ + 2P«)P - p6 - 2P^ 



(120) 



33 




FIG. 21: Susceptibility hexagons on either H or B lattices. 



Now we consider the diagrams of Fig. |2T1 We have 



AR{a) = 2yi^Rl_.^[P\l + 2R + 3R^) +4PR^ + K^] 

+2Rk-^R[P(l - P - R){l + 2R + m^)] + Rk-^R^ 
+R^[{1 -P- Rf{l + 2R + 3i?2) + 4(1 - P - R)R^] 



(121) 



AR{b) = -AyR'[R{l - P) + Rk-sP] 



(122) 



AR{c) = -4yi^RRk^3[{l + R + R'^ + R^ + R^)P] + R^[l - P(l + R + R^ + R^ + RX^) 



AR{d) = -Ay 



Ril -P) + Rk-^P 



R[l - P{l + R + R^ + R^)] + Rk-3Pil + R + R:' + R-') 



2 I d3\ 



AR{e) = -Ay 



R[l - P(l + R)] + Rk-sPil + R) 



R[l -P{l + R + R^)] + Rk-sPil + R + R' 

AR{f) = 2yR' 
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(124) 



(125) 
(126) 



(a) 



-2 




(c) 







FIG. 22: Susceptibility hexagons on either H or B lattices. 
AR{g) = 4yi^RRk-3[{l + R + R^ + R^)P] + R^[l - P{1 + R + R^ + (127) 



AR{h) = Ay 



R{1 -P) + Rk^sP 



R[l - P(l + i? + R^)] + Rk-3P{1 + R + R^ 



(128) 



R[l - P(l + R)] + Rk--iP{l + R) 



(129) 



We call the sum of all these contributions Vh b- We have 



Th,b = 4a V Rl_, 4R^[P - P^] + i?M + Rk- 



R*{P'^ - P) - 4P^P + R^ 



AR^P^ + AR^P + 4PP6 - 2PS . 



(130) 



Next we consider the diagrams of Fig. |221 

AR{a) = 2y|p^_3[P^ + 4P^P + 2p2p2^i?p2j 

+PPfc_3[2PP(l - P - P) + 4p2p(l -P-R)+R^] 



(131) 



AP(6) = -4?/P^ 



P[l - P(l + P)] + Pfc_3^(l + R) 



(132) 
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FIG. 23: Susceptibility hexagons on either H or B lattices. 



Ai?(c) = -2yR'[Ril-P) + R,_sPf 



(133) 



AR{d) = -AyR^ 



R[l -P{l + R + R^ + R")] + Rk-zP{l + R + R' + R')] 



2 , d3^ 



(134) 



Ai?(e) = -AyR 



R{1 -P) + Rk-^P 



R[l -P{l + R + R^)] + Rk-^P{1 + R + R^)\.2b) 



AR{f) = -2yR 



1 2 



R[l - P(l + R)] + Rk-sPil + R) 



(136) 



AR{g) = AyR' 



R[l - P(l + + + Rf,_^P(l + R + R^) 



(137) 



AR{h) = AyR 



R[l - P{1 + R)] + i?fc-3P(l + R) 



R{1 -P) + Rk-sP 



(138) 



AR{i) = 4yR^[R{l - P) + Rk-sP] 



(139) 



We call the sum of these contributions Thc- We find that Thc = b- 



Next we consider the diagrams of Fig. |2S1 



AR{a) = ?/|i?^_3[i?^ + AR-^P + 2R'P'] + RRk-^[AR'P{l - P - R) + R^ 
+p2[2(l - p - p)2p2 + 4(1 _ p _ p)p3] I 



(140) 



AR{b) = -AyR^ 



R[l - P(l + + R")] + Rk^-sPil + R + R^ 



(141) 
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AR{c) = -AyR^ 



R[l - P(l + R)] + i?fc_3P(l + R) 



[R{1 - P) + Rk-3P] (142) 



AR{d) = 2yR^[R{l - P) + P,„3P] 



(143) 



AR{e) = AyR^ 



R[l-P{l + R)] + Rk-zP{l + R) 



(144) 



We call the sum of these contributions Thd. We have Th d = Vh^b/'^- Then the total 
contribution to the renormalization of R from hexagons is 



APnex = 2p^aHRk-i 



5R\P - + r5 



20P^(P - P2) + 



- 25P^(P2 - P) 



+P, 



k-3 



50(p2 - P)P^ - 22PP^ + 



+ 22P^P-11P^-2P^ 



(145) 



VII. CONCLUSION 



We now have all that we need to evaluate whether or not the perturbed ratio r, i.e. 
the right-hand side of Eq. (fTIj) . is unity or not to order Using Eqs. ()145|) and (fTOj) 

we see that the contributions from the hexagons do not contribute to the rhs of Eq. (j74j) . 
Furthermore, using Mathematica Js^, we summed the contributions from the squares in 
Eqs. (03), (EH), dlHI), and dTTTHl . which led to the resuh (at the transition in the EOS) that 



1 + (3ap72)p2[i? _ Pfc_3][(2 - k)R + aQRk--i] + Oia 



(146) 



In writing this result we have assumed that the effect of l/rf perturbations is merely to 
shift the location of the pole in the susceptibility. It is not obvious that this is a correct 
assumption, since we do not know a priori that for finite d the suceptibility diverges (as 
it does in the other problems for which we cited the use of the 1/d expansion) [31.]. A 
discussion is given in Appendix B to justify this assumption. 

Since p ~ l/o", the term proportional to is of order a~^. Note that to leading order in 
1/(T we have that 



R„ 



(147) 



which implies that (2 — k)R + aQRk-3 = 0{l/a), so that to order (1/cr)^, the result of 
Eq. ()146|) is that r = 1, which means that to this order the singularity in the EOS and the 
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divergence in the susceptibility coincide. Note that this result is not a trivial one in that 
it implicitly involves Eq. ()147|) . Also, it can hardly be a coincidence that the subtleties of 
diagram counting lead to this result. It would not totally surprise us if this result could be 
obtained to all orders in 1/a by some type of Ward identity, even in the absence of a field 
theoretical formulation. 

The coincidence of the two singularities suggests that the unusual nature of the k > 3-core 
transition survives in large, but finite, spatial dimensions. Given the absence of numerical 
confirmation of such a transition in simple isotropic fc-core models on square, cubic, and 
triangular lattices, our results motivate further numerical studies of models that are more 
mean field-like in the sense that the range (number) of nearest neighbors a + 1 is larger 
than previous models. Of course, to compare with our results one would need to study 
hypercubic lattices at large d. Since "culling" leads to nonlocal truncations, it will not be 
easy to simulate systems at large d to check our work. However, should such a program be 
undertaken, one might want to have actual expressions for the shift in the critical point and 
the shift in the jump of the order parameter due to finite dimensionality. To get the shift 
in one needs to evaluate the rhs of Eq. ()24|1 . using Eqs. ()43|1 and ()7()|1 . This result will 
include all corrections up to and including order However, it is necessary to expand 

the quantities such as po, Qo, and Rm in powers of 1/cr, as is indicated by Eqs. (fT7|) and 
fl29|) . Furthermore, to get the shift in the jump in the order parameter involves evaluating 
the rhs of Eq. (j26|) and for this we had recourse to Mathematica js^ . 

In order to assess the implications of our work when placed in context with known results 
we refer to the phase diagram in the k-d plane in Fig. |^ where we indicate existing results 
for k-core percolation. The two regimes which are most securely established are a) for = 2 
the model has js^ the same critical point as ordinary percolation and b) for k > d + 1 the 
critical concentration for percolation has been shown 3| to be Pc = 1- The present work 
suggests that the hybrid transition found for d = oo persists into the regime k <^ d < oo. In 
addition, numerical work indicates that the transition may be continuous for k = d = 3 Q| . 
Since this is the only data point which gives a continuous transition for A; > 2, it is important 
to confirm this result on larger samples, although our discussion will assume the validity of 
this result. The 1/d result of the present work indicates that the hybrid transition occurs 
for A; > 3 at large d, so it is clear that the continuous transition seen ai k = d = 3 
must disappear as d is increased for k = 3. The details of this boundary between hybrid 
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and continuous transitions is unclear at present, except that it seems almost certain that 
/c = 3 is essentially different from k = 2. Furthermore, at present there is no evidence of 
yet another phase between the hybrid phase and the Pc = 1 phase at high d. This phase 
boundary might remain at = + 1. Alternatively, since finite-dimensional fluctuations 
reduce Pc in the hybrid phase from its mean-field or Bethe lattice value, it is possible that 
the phase boundary between these two phase falls below the line k = d+ 1, as indicated by 
the question mark in Fig. |211 From the perspective of Fig. |^ it would be interesting to 
develop a realization of k-core percolation for noninteger k. 

It is interesting to note how the l/cr ex pan sion for fc-core percolation compares to that 
for other systems. For self-avoiding walks j26|, the Ising mo del H, 

spin glasses and 



for percolation |29|, the expansions for the critical value of the coupling constant involves 
coefficients of l/cr", which are rational fractions (i. e. the ratio of two finite integers). In 
contrast, for lattice animals one sees the appearance of the transcendental number e in the 



coefficients of l/cr"' j2^. From Eq. (20) one sees a similar result for fc-core percolation. It 

n 

has been noted 1171 that the unusual mean-field value of the correlation length exponent for 

n 

k-coie percolation is identical to that for lattice animals j2^. Thus, it is not surprising that 
these two models show similar unusual characteristics in their 1/a expansion for the critical 
coupling constant. 



VIII. DISCUSSION 

The model of fc-core percolation sets a constraint on the number of occupied neighboring 
bonds (or sites in the corresponding site problem). So the physics of fc-core pertains to 
systems with nontrivial constraints such as the onset of long-range orientational order in 
solid ortho-para H2 mixtures and the onset of rigidity in a mechanical system. Given the 
equivalence between rigidity percolation and {g + l)-core percolation on the Bethe lattice 



12| . the results here suggest that provided g is at least 2 in high dimensions, the percola- 
tion transition should be of a hybrid nature where there is a jump in the order parameter 
accompanied with a diverging length scale. We note that there exists a field theory of 
rigidity percolation that exhibits a purely first-order transition, which is very different from 
the hybrid nature of k-core percolation. Given our results, as well as the results of SLC, it 
seems that this task of constructing a field theory for rigidity percolation should be revisited. 
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FIG. 24: Summary of known results for fc-core percolation on the d-dimensional hypercubic lattice. 
Triangles represent ordinary percolation for k = 1 and biconnectedness for k = 2 (whose critical 



exponents are governed by the fixed point of ordinary percolation |3^). Circles represent systems 
for which pc was shown to be unity for d = 2 by van Enter and by Toninelli et al. 19] for 
larger d. The filled square represents the result of Ref. for d = k = ^ where the transition 
appears to be continuous. This result therefore suggests the existence of a region labeled "Cont" of 
a continuous transition with a critical point probably similar to that of ordinary percolation. The 
filled diamonds represent the present work suggesting that the hybrid transition survives for large 
d (but with k not comparable to d). The question marks indicate that the boundaries between the 
continuous transition (filled square) and the hybrid transition, as well as that between = 1 and 
the hybrid regime, are unclear at present. 



le^dynamics of an interface separating 
When the disorder is strong, it is 



Other systems with nontrivial constraints include t 
two magnetic domains in a random field Ising model 
likely that more than several neighbors must have flipped previously in order to propagate 
the interface, again, giving rise to /c-core physics. However, finite-dimensional simulations of 
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such interfaces have only uncovered ordinary percolation exponents so far [36[, as opposed 
to some sort of hybrid transition. 

While not directly related to fc-core percolation, spin glasses are systems where a set of 
nontrivial constraints cannot all be satisfied. This property is otherwise known as frustration. 
Curiously enough, it turns out that the p > 3-spin glass model also exhibits the unusual 
hybrid transition in mean field with the same exponents as fc-core |37]. It would also be 
interesting to see if a l/ci expansion for this system yields the same results that we find for 
A;-core. 

The physics of glass-forming liquids, a.k.a. the glass transition, has received a lot of 
intense investigation over the years yielding a wide range of approaches js^- One approach 
to modelling glassy dynamics is based on kinetically constrained spin lattice models, like 
the Frederickson- Andersen (FA) model js^. In this model, down spins on a lattice denote 
coarse-grained regions of high mobility of the liquid, while up spins denote regions of low 
mobility. A magnetic field favors up spins creating large regions of low mobility. As the 
temperature of the system is lowered, more and more regions of low mobility are created 
eventually leading to kinetic arrest of the liquid over large time scales. An important kinetic 
constraint on the motion of the spins is introduced where a randomly selected spin can flip 
only if the number of neighboring downs spins is equal to or greater than some integer whose 
maximum is cr + 1. This constraint models the caging effect observed in glass- forming liquids 
where particles become trapped by transient cages made up of their neighbors. 

There is an exact mapping between fc-core percolation and the clustering of low mobility 
regions in the FA model (40, |4l|. The mapping to fc-core is not unexpected given that 
up/down spins can be mapped to occupied/ unoccupied sites, and therefore, the kinetic 
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for details. Given this mapping. 



constraint maps to the fc-core condition. See Refs. 
our conclusions apply to that type of lattice model, and in high d one should observe a hybrid 
transition at finite temperature, provided k is not comparable to d. We note that this regime 
may not be directly relevant for the glass transition due to the existence of finite clusters, 
i.e. the cube in d = k = 3 j^]. Then the density of particles in the liquid phase is not 
homogeneous. However, recent experiments on metallic glasses show short-range icosahedral 
order even in the hquid phase so it is not clear if the density should be assumed to be 
homogeneous, at least over some time scale 4^. Therefore, a thorough understanding of k- 
core beyond mean field is of the utmost importance to understanding such transitions as the 
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FIG. 25: Enumeration of the susceptibility diagrams for the Bethe lattice. Bonds which are parallel 
involve the same local displacement and sites which are close to one another coincide when the 
diagram is implemented on the hypercubic lattice. 

glass transition and the onset of orientational ordering in solid ortho-para H2 mixtures, and 
our work seems to be the first analytic result to deal with its finite-dimensional fluctuations. 
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APPENDIX A: COUNTING SQUARES 

Here we consider in detail the counting of squares with excluded volume corrections which 
represent corrections of order 1/a relative to leading corrections due to insertion of squares. 
More specifically, we consider the enumeration of pseudo square insertions at a single vertex, 
like those shown in Fig. |2H1 Here there are two distinct types of square insertions, all of 
which have to subtracted off because their topology is different on a hypercubic lattice than 
on a Bethe lattice. In panel (a) the situation is simple, in that no sides of the square are 
equivalent to the other bonds even on a hypercubic lattice. The diagram in panel (b) has 
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to be subtracted off, but note that because two bonds coincide on the hypercubic lattice, 
this diagram is in the class of diagrams of the topology of Fig. ^[ which here is shown in 
panel (c). To avoid subtracting this diagram twice, we do not allow it in the enumeration of 
diagrams we wish to subtract off here. Thus we wish to enumerate pseudo squares (diagrams 
which would involve a square on the hypercubic lattice), but only those in which all bonds 
are distinct from the bonds of the chain into which the square is inserted. If the two bonds 
{u — S and S — v) entering the vertex are in the same direction there are of order ways 
to form the square. So the number of configurations of this type is of order a^. We next 
count the number of configurations when the bond S — v is perpendicular to the incoming 
bond. There are then a — 1 ways to choose this bond. As shown in panels (d), (e), and (f), 
the square can then either (d) be in the same plane as bonds u — S and S — v (and there 
are two ways to do this), or (e) the square can have two bonds parallel to either bond u — S 
or bond S — v (there are four ways to do this) and then there are cr — 3 ways to choose the 
other bonds of the square to be perpendicular to bonds u — S and S — v, or (f ) all bonds of 
the square can be perpendicular to the bonds u — S and S — v (there are (a — 3) (a — 5) way 
to do this. So in all, the number of configurations, Ng, of the square and of bond S — v is 

Ns = [a^ + Oia)] + (a - 1)[2 + 4((t - 3) + (a - 3)(a - 5)] 

= a'[l + {A/a)] + 0{a) . (Al) 

This justifies the prefactors for diagrams (b) and (f) of Fig. ^| 

The analysis for the diagrams of Fig. El is similar. In panels (a) and (b) of Fig. 1^ we 
show diagrams which have a different topology on the Bethe lattice than on the hypercubic 
lattice and therefore which must be subtracted off. In panel (c) we show a special case when 
bond Q-R coincides with a bond (U-w) in the chain. As before, although we do need to 
subtract off this diagram, we do not include it in the present enumeration because it is of the 
topology of panel (d) which is included in our analysis of Fig. If bond S* — T of panel a 
is parallel to the incoming bond then there are cr — 1 ways to complete the square and a — 1 
ways to add the bond exiting the square. When bond 5 — T of panel b is perpendicular to 
the incoming bond (there are a — 1 ways to do this), there are a — 2 ways to complete the 
square and then a — 1 ways to add the bond exiting the square. Thus in all, Ns, the number 
of ways of configuring the square and the exit bond, U — w, is 

= (a-l)2 + (a-l)((T-2)(a-l) 
43 
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FIG. 26: Enumeration of the susceptibility diagrams for the Bethe lattice. Bonds which are parallel 
involve the same local displacement and sites which are close to one another coincide when the 
diagram is implemented on the hypercubic lattice. 

= a^[l-{'i/a)] + 0{a) . (A2) 

This calculation justifies the prefactors for diagrams (b) and (e) of Fig. El 

The diagrams of Fig. ^Jare more straightforward because for them one does not have to 
consider the possibility of double counting the subtractions. 



APPENDIX B: SHIFT OF SINGULARITY IN THE SUSCEPTIBILITY 

In Sec. VI we evaluated the ratio r in the series for the susceptibility and interpreted the 
result as giving a shift in the singularity. Here we justify this interpretation. Note that the 
vertex renormalization (in which the unperturbed ratio tq is replaced by tq + tqA) gives the 
susceptibility as 

X = E^o[l + ^A]. (Bl) 

n 

In writing this result we noted that for a chain of n bonds the vertex renormalization could 
be placed at any one of order n vertices. Of course, there are end effects, so that really, if 
we include end effects, we would write 

X = E^o[l + ^A + e] , (B2) 

n 
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where both A and e include all contributions of order o""^ and a~^. Now we consider the 
contributions of two vertex renormalizations. If the two vertices do not interfere with one 
another then their contribution is tqA^. For a chain of n bonds, there are of order such 
configurations. So in analogy with Eq. ()B2|) . these contributions are of the form 

5X = E^oKAV2 + ne + r/] , (B3) 

n 

where e comes from configurations involving two squares or higher order configurations, none 
of which were counted up to order a^^ in Eq. ()B2j) . These contributions involve either one 
square (or hexagon) near an end point and the other in the interior of the chain, or two 
interfering structures in the interior of the chain. Also 1] comes from configurations in which 
all insertions are near an end of the chain. 

It seems reasonable to assume that the dominant contribution to x can therefore be 
written as 

X = E ^0 [1 + + AV2 + ...]- E ^^ (B4) 

n 

where r = ro[l + A], as we found in Sec. VI. The analysis of this appendix indicates that is 
appropriate to identify r as the renormalized ratio of a geometric series. 

One might ask whether this identification is unique. Could the result of Eq. ()B4|) arise 
from a rounded transition for which 

X = ^E^o[(l + «r + (! + «*)"] (B5) 

^ n 

for a suitably chosen value of the complex-valued parameter a? This form of x yields 
(keeping only relevant terms) 

X = E ^0 [1 + + «*)/2 + ln\a' + a*')] . (B6) 
For this to be of the form of Eq. ()B4|1 we find that a must satisfy 

\ia' + a*') = ^[a + a*]\ (B7) 

which implies that a* = a. Therefore, we cannot have a rounded transition by having the 
real valued critical point for the Bethe lattice split into a complex conjugate pair of critical 
points (which would give a Lorentzian susceptibility with a width of order o"~^). So the 
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form of Eq. ()B4|1 is only consistent with a shifted pole in the susceptibility, as we implicitly 
assumed in Sec. VI. 
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